Bulletin of Electrical Engineering and Informatics 
Vol. 13, No. 2, April 2024, pp. 893~904 
ISSN: 2302-9285, DOI: 10.1159 1/eei.v13i12.5815 o 893 


Simplifying the electronic wedge brake system model through 
model order reduction techniques 


Mohd Hanif Che Hasan!, Mohd Khair Hassan?, Fauzi Ahmad’, Mohammad Hamiruce Marhaban?, 


Sharil Izwan Haris‘, Ehsan Arasteh? 
'Department of Electrical Engineering Technology, Faculty of Electrical Technology and Engineering, Universiti Teknikal Malaysia 


Melaka, Melaka, Malaysia 


Department of Electrical and Electronic Engineering, Faculty of Engineering, Universiti Putra Malaysia, Serdang, Malaysia 
3Department of Automotive, Faculty of Mechanical Technology and Engineering, Universiti Teknikal Malaysia Melaka, Melaka, 


Malaysia 


Kolej Kemahiran Tinggi MARA Masjid Tanah, Melaka, Malaysia 


‘Department of Mechanical and Aerospace Engineering, University of California, California, United States 


Article Info 


ABSTRACT 


Article history: 


Received Feb 2, 2023 
Revised Sep 12, 2023 
Accepted Oct 3, 2023 


Keywords: 


Brake by wire 
Electronic wedge brake 


The electronic wedge brake (EWB) uses self-reinforcement principles to 
optimise stopping power, but its mathematical model has various actuation 
angles and system dynamics making controller design complex and 
computationally burdensome. Therefore, the model order reduction (MOR) 
is made based on three factors that may have a negligible influence on the 
EWB system: the motor inductance, lead screw axial damping, and wedge 
mass. Six reduced order model (ROM) types were proposed when one, two, 
or all factors were ignored. The ROM accuracy was analysed using the 
frequency and time domain. The percentage of root means square error 
(RMSE) response value between the EWB benchmark model, and the 
predicted response based on the ROM was found to be less than 2%, with 


Model order reduction ROM size reduced from 5 to 2 orders. It guarantees that the new ROM series 

Modelling will be useful for simpler EWB controller design. The proposed ROM 

Reduced order model simplifies the original model drastically while retaining accuracy at an 
adequate level. Even though the simplest EWB model is a 2™ order linear 
system, the best ROM vary depending on EWB design parameters. 
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1. INTRODUCTION 

The brake by wire (BBW) is a next-generation brake system that entirely replaces the traditional 
brake mechanism. The electro-mechanical brake (EMB) is the most accurate description of BBW because 
electronically powered actuators generate the braking force. The EMB has the following advantages: it is 
noiseless, has a faster response time, is easy to create a control system, and is more ecologically friendly [1]. 
Most EMBs have an electric motor, a reduction gear, a floating disc brake calliper, and a 24/42 V power 
supply [2], [3], and consume a lot of energy [4]. Since existing automobiles only use a 12 V power supply, 
the electronic wedge brake (EWB) technology has been a popular alternative because it can run at the same 
voltage as the power source and uses around one-tenth of the power [5]. The EWB is from the EMB category 
that, by employing the self-reinforcement principle, can reduce energy usage. 
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Two types of EWB designs are categorised based on their actuation angle. The first type is the 
tangential angle, where the actuation angle is parallel to the disk surface. For this configuration, the actuation 
angle is 0 degrees. Therefore, the clamping force acts perpendicularly to the actuation force produced by the 
motor. In comparison, the second type is the optimised angle, when the actuation angle is set at the same 
angle as the wedge angle [6]. As discussed in [7]-[11], different actuation angles are analysed using 
analytical and experimental approaches. It is found that the best driving angle is when it is set equivalent to 
the wedge angle. 

The EWB system can be modelled using physical modelling or system identification (SI). The SI 
method has been used to develop dynamic models of EWB such as Box-Behnken [12]. The physical 
modelling was first used in [6], where a state-space representation of the EWB model was used to construct a 
braking force controller. The EWB model was mathematically derived using ordinary differential equations 
(ODE) and was developed in three separate sub-assembly parts: the DC motor, lead screw, and EWB brake 
heart. As the investigation progresses, subsequent numbers of mathematical modelling of the EWB system 
have been developed [13]-[20]. In the EWB model, five states are identified: motor current, motor angular 
velocity, angular motor position, wedge speed, and wedge position. Three of these state vectors are derived 
from the motors and lead screw dynamics, while the other two are from the wedge model. 

Furthermore, Hwang and Choi [21] describes the EWB model for tangential actuation type from an 
electric motor to a wedge mechanism in detail, with all differential equations for all components combined. 
However, the derived model was intended for basic EWB, which uses an actuation angle parallel to the disk 
surface (tangential type). Although the second-order system was widely used to describe wedge mechanism 
dynamics, there was also an 8" order system introduced to produce a model with at least 10% higher 
accuracy than others [22]. For this consideration, the entire EWB model mathematical description has up to 
the 13" order system after combining an electric motor (3™ order) and a roller screw (2™! order) dynamic. 
However, researchers did not use this model because of its enormous size and many parameters. 

Due to the different types of actuation designs will have different model structures, the first flexible 
version of the EWB model is introduced, which can be used for two famous EWB designs based on 
tangential and optimised actuation angles [23]. The new generalised model is more compact and simpler than 
the various versions of EWB while retaining model accuracy. Despite the 5" order linear model being 
capable of describing the EWB system dynamic well, some nonlinearity is neglected in this model, such as 
gear backlashes, coulomb frictions, and disk gap clearance inherited from conventional EMB [24]. 

Model-based approaches, which are prevalent in modern control strategies, benefit significantly 
from representation in state space. Modern control strategies applied on EWB so far such as proportional- 
integral-derivative (PID) [2], [6], linear quadratic regulator (LQR) [25], Youla parameterisation [26], sliding 
mode control (SMC) [17], [21], and active disturbance rejection control (ADRC) [17] is highly dependent on 
a model that represents the actual model. The model order reduction (MOR) simplifies controller design for 
complex models. It falls into two categories: conceptual and mathematical/numerical/data-driven [27]. Model 
reduction of linear-nonparametric dynamical systems has reached a considerable level of maturity. However, 
parametric model reduction has emerged more recently as an important and vibrant research area [28]. Even 
though research in parametric model reduction now mainly concentrates on nonlinear, parameterised, and 
coupled problems, linear problems are still some unexplored challenges, especially in certain applications. 

Reducing the complexity of the 5 separate states involved in designing a controller for EWB is 
difficult due to limited sensor access. However, neglecting certain factors such as motor inductance, lead 
screw axial damping, and wedge mass can simplify the 5" order linear EWB model to a 2 order linear 
system [21]. Nonetheless, no further research is done to see what happens if not all these variables are 
considered. Based on these three factors, we performed 2 cases based on our EWB parameters and realised 
that the number of factors that can be ignored depends on EWB design parameters [29]. Motivated by this 
result, we have expanded our study to several cases, as in this article and provided detailed derivation and 
assumptions. We proposed 6 cases with different EWB simplified models. Based on our EWB design, further 
analysis is made to evaluate each model accuracy, with the flexible, 5 order linear model being used as a 
benchmark. All the symbols and parameters used in this paper are based on previous work reported in [23]. 

There are four sections in this article. The first segment focuses on the EWB modelling review of 
other relevant work. The equations involved in complete parametric modelling of the EWB are then listed, 
followed by the 6 cases based on different scenarios. The simplified model for each case is proposed, and the 
output response of the model is simulated and observed based on step response (time-domain) and bode plot 
(frequency-domain) form. Besides, the root means square error (RMSE) value is used to find the best- 
simplified model before ending with the conclusion in the last section. 
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2. METHOD 
2.1. Original model 

As shown in Figure 1, the brake actuator consisted of the DC motor and the roller screw. The DC 
motor actuator was used in this system to move the input torque and force the centre brake through a roller 
screw, which transforms the angular motion of the DC motor into the axial motion of the coil in the centre of 
the brake. On the other hand, the heart of the brake, consisting of a wedge mechanism, a calliper, and a brake 
pad, may have diagonal motions, providing a clamping force to the brake disk. 


y Roller ee | p 


Figure 1. EWB schematic diagram 
There is a total of 7 equations that represent the dynamic of the EWB system. Several terms are 
introduced to facilitate the derivative process as follows: 
a, = Koa tana (tan a — u) 


az = My (tan? a + 1) 


LaNa 
a3 = 21 
a, = cos fp 
Where 
f= , for Tangential Actuation EWB 


æa for Optimised Actuation EWB 


The complete model of EWB consists of a few physical parameters. The six parameters of DC 
motor are resistance (Rm), inductance (Lm), electromotive force constant (Ke), torque constant (Ki), inertia 
(Jm), and viscous friction constant (Dm). The lead screw consists of 5 parameters which are reduction gear 
ratio (Na), steadiness (Ka), viscous damping (D,), pitch (La), and efficiency (n). The wedge mechanism and 
brake pad parameters are wedge mass (Mw), wedge angle (a), actuation angle (P), calliper stiffness (Kca), and 
brake pad coefficient (u). The DC motor electrical and mechanical dynamics describe by (1) where Im is a 
motor current, Om is a motor angle, and œm is an angular speed of a DC motor. 


vee Me. a = 
a i Wm = Im + 7 Vin (1) 
5 D. K, 1 
Wm = -7 Om + Tlm = 7, [screw (2) 
m m m 


For the lead screw dynamic that connects DC motor and wedge, the motor force (Fm) and screw 
torque (Tscrew) can be defined as (3) and (4): 
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Xw * Vw 
Fn = Ka (250m, — 72) + Da (236, - *) (3) 
Fm 
Tscrew = a3 n (4) 


Where X, is a displacement and V,, is a speed of a wedge. By substituting (3) into (4), the torque screw 
(Tscrew) can be defined as (5): 


Xw - W 
T screw = = [Ka (436 cs že) + Da (436), = “| (5) 
For wedge mechanism dynamic is described as (6): 


Ve = (-2) Xu + (=) Fn (6) 


a2 


Note that these six equations used are the simplified version when the angle of actuation (B) can be either 0 or 
equal to wedge angle (a) as in [23]. Finally, the clamping force (F.) of the wedge depends on the wedge 
displacement (Xw) in the x-direction, calliper stiffness (Kca), and wedge angle (a). 


F = KeaXw tan a (7) 
2.2. Assumptions 

The first assumption that can be made is the motor inductance factor. Motor inductance usually has 
a minimal value, making it challenge to be a dominant pole. When this factor is ignored, (1) can be rewritten 
as (8): 


Ke 1 
Im = = p Om + 2 m (8) 


By substituting (8) into (2), both of DC motor equations is combined and become a single description. 


Wy = — (FE SEE) tom (T Un — > Trove (9) 


Jm  JmPm JmRm 


Meanwhile, the second possible assumption is based on the axial damping factor. In most cases 
where the mechanical component is so stiff, the axial damping only has little effect on the system response. 
Therefore, the equation for the axial connection in (3) and (5) can be rewritten as (10) and (11): 


En = (Kas) Om — ($) Xw (10) 
Tocrew = (E) Om - (2) Xy a1) 


The third assumption, on the other hand, is based on the mass of the wedge. In a typical EWB 
mechanical system, the light mass wedge is usually coupled with a high degree of stiffness. Thus, resulting to 
a stiff mechanical part with very high bandwidth. So, the performance of the system response largely depends 
on the electrical actuation system. Accordingly, the dynamic of a wedge is ignored by setting the mass of the 
wedge to zero. So that (6) is transformed into: 


Fn = (a,a4)X,, (12) 

In this study, 6 cases were selected for further evaluation. These cases were derived based on 
assumptions and resulted in state-space representations. The first 3 cases are based on a single factor being 
ignored: motor induction, axial damping, or wedge mass. The fourth and fifth cases are combinations of two 


factors being neglected. The sixth case considers a combination of all factors being neglected. 


2.3. Case 1—neglection of DC motor inductance 
In (13) yields when (5) is substituted into (9). Meanwhile, substituting (3) into (6) yields (14): 
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- — (Eam )x + (2em yy - (=) g — (2z + KtKe + Pass) + ( Kt yy 13 
Om fren y: NJma4 w nJm i Jm JmRm nJm Om JmRm m, ( ) 


Wy a a a ae) mt Ce) oe 


Using (7), (13), and (14) will reduce the order of state-space and hence, yields: 


x = Ax + Bu 
y = Cx + Du (15) 
x = [Xw Ww. 0m, wml” 
yak 
0 1 0 0 
_ Kata1a4? _ Da Kqa3 Daa3 
A= a2a42 az2a42 a2a4 a2a4 
0 0 0 1 
Kaas Daas _ Kaas? _ Dm _ Keke _ Dats? 
NJma4 NJma4 nJm Jm  IJmRm nJm 
T 
: Koai tan(a) 
B=| o |, andC= 0 
Kt 0 
0 


JmRm 


Therefore, system order reduces by 1 when neglecting motor inductance. The EWB model then 
becomes a fourth order system with four state vectors instead of the fifth order in its original model. 


2.4. Case 2-neglection of axial damping 


In the 2" case, (16) yields by substituting (10) into (6). Furthermore, substituting (11) into (2) 
produces (17). 


Vy = — (Eie) x, + (E2) Op (16) 
tin = (EE) Ky — (EE) Bon — (JE) om + (FE) In a7) 


Finally, by combining (1), (7), (16), and (17), the complete model of the EWB actuator without axial 
damping dynamic factor is obtained: 


x= [Xw Vive 0m, Ons Lal? (18) 
0 1 0 0 0 
Ka+a1a4? Kaa3 
T az2a42 0 a2Q4 0 0 
A= 0 0 0 1 0 
Kqa3 0 _ Kaa? _ Dm Kt 
NJma4 nJm Jm Jm 
0 0 0 2e aft 
Lan La 
A Kear tan(a) i 
0 
B= 0 , andC= 0 
y 0 
La 0 


In this case, even though the system is simplified, the system order remains like the original model. 
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2.5. Case 3-neglection of wedge mass 
Negligence for case 3 is possible when the wedge mass is very light. From (12) is substituted into 
(4) and then substituted into the motor in (2) will yields (19). Meanwhile, (12) with (3) will produces (20). 


. (419304) y _ Dm Ke 

Wm = ( Los F int ha (19) 
3 2 

X, = — (FEE) xy + (AE) On + (stg) (20) 


By combining (1), (7), (19), and (20), the reduced order model (ROM) for case 3 is produced. 


x = [Xw 9m Wm Im] (21) 
Ka+ 2 Ke 
— a n azđ4 0 
0 0 1 0 
A= 
= fsa. 0 _ Pm Ke 
nJm Jm Jm 
0 0 nie, yee 
Lin Li 
: Kca tan(a) a 
B=|Qg|, andC = 0 
i 0 
Tie 0 


Therefore, system order reduces by one when considering wedge mass is very light. The EWB 
model becomes a fourth-order system with four state vectors instead of the fifth order in its original model. 


2.6. Case 4—neglection of motor inductance and axial damping 


In case 4, (11) is substituted into a simplified equation for neglecting motor inductance case as in 
(9), which then produces (22): 


tin = (pes) Hw (Cia) Om Ga t Tana) om + GE) Mo 22) 


Using (7), (16), and (22), the ROM size of the state-space yields (23): 


x= [Xw Vv, 9m, ml” (23) 
0 1 0 0 
Ha Ka + nae 0 Kaa3 0 
a2A,4 a2A,4 
hee 
0 0 0 1 
Kaz 0 K,a3” Dm K,Ke 
Mm, Wm Jm JmRm 
0 T 
0 Kca tan(x) 
B=| 0 |, andC= 0 
K; 0 
0 
ImRm 


2.7. Case 5—-neglection of axial damping and wedge mass 
For the 5" case where axial damping and wedge mass are ignored. In (10) with (12) produces (24). 
Meanwhile, (25) yields by substituting (22) into (17). 
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Xy = [pete] 6, (24) 
tin = (Fass T ma) Pn Gn) me + Ga) 5) 


Substituting (24) into (7) generates: 


F, = Koai tan a ($) Om (26) 


Ka+a1a4? 


Finally, by combining (1), (25) and (26), the complete model of the EWB actuator when axial damping and 
wedge mass are neglected is obtained. 


x= [Om wm, Im] (27) 
0 1 0 
Ka?a3?  _ Kaag? _ Dm Kt 
A= NJm(Ka+a1a4?) nJm Jm Jm 
0 Ke Rm 
Lah Lm 
0 Keqi tan(«)Kga3a4 T 
_10 = Ka+a1a4? 
B= AE and C = 0 
Lm 0 


Therefore, in this case, the system order reduces by 2. The simplified model of EWB becomes a 
third-order linear system. In this case, an estimated transient behaviour purely comes from motor behaviour. 


2.8. Case 6—-neglection of all 3 factors 


The maximum reduction of the model can be made by assuming that all of the three aforementioned 
factors can be ignored. Substituting (24) into (22) produces the single as (28): 


Bo tae Kqa3 e aag? z)|0 ae + Ake) ( Kt \y, 2 
Om [fe Kqta,a4? us Jm  JmPm Oma JImRm/ ™ (28) 


Using (26) and (28), the EWB ROM in the state-space is obtained. 


x = [Om Om)" (29) 


0 1 
A= S$ Kaaz? ( a1a4? ) _ Dm AS KtKe 
nJm \Kat+a1a4? 


0 K cal tan(X)Kga3a4 T 
B= | Kt | and C = | Ka+a1a4? 
0 


Jm Jmm 


ImRm 


Thus, the system can be represented only with a motor angle (Om) and motor angular speed (@m) in 
the simplest model. The system order is reduced from five to two. 


3. RESULTS AND DISCUSSION 

Based on derived state-space via MOR discussed before, it concludes that the size and complexity of 
an EWB model are possibly be reduced by using 3 possible assumptions. The original model which is a 5" 
order system can be reduced as low as a 2™ order system and summarised in Table 1. The best ROM is 
chosen for its accuracy and simplicity. Assumptions for controller design are determined. Frequency and 
time-based analyses were performed using a bode plot and step response. 5 order EWB model was used to 
evaluate the accuracy of simplified models. 

Frequency-based analysis of six ROMs is plotted in two bode plot forms as in Figures 2 and 3. From 
both figures, only several cases are considered acceptable. Only cases 2, 3, and 5 were permitted while 
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showing results close to the benchmark model. The case 2 bode plot follows the benchmark model closely 
when compared to others for both the magnitude and phase values. Cases 3 and 5 also follow the benchmark 
model closely at low-level frequency. However, both models start showing different behaviours at a high 
frequency that starts from 30 k rad/s and above. 


Table 1. The EWB ROM in state space summary 


Case no. Neglecting factor Independence states _ Clamping force (F,) 
1 Motor inductance (Lm) Xw Vr Om, Om Obtain from Xy 
2 Axial damping (D,) Xw Vr Om: Ome Im Obtain from X,, 
3 Wedge mass (Mọ) Xw Om Wm Im Obtain from X,, 
4 Motor inductance (Lm) and axial damping (Da) Xw Vy, Om: Om Obtain from X,, 
5 Axial damping (Da) and wedge mass (My ) Om Wm» Im Estimate from 0, 
6 Motor inductance (Lm), axial damping (D,), and wedge mass (My) Om, ®m Estimate from 0, 


Bode Diagram 
100 T T T T T T 


-100 


EWB 5" Order Model 
— — -Case 1: Lm Neglected 
I< XCase 2: Da Neglected 
pee kee Case 3: Mw Neglected 


Magnitude (dB) 


-200 F 


Phase (deg) 


-450 & pil n soul l i METTETE ENET 


1071 10° 101 10? 10° 104 10° 10° 
Frequency (rad/s) 


Figure 2. The bode plot of the EWB original model versus first 3 cases of ROMs 


Bode Diagram 


-100 


EWB 5!" Order Model 
— — -Case 4: Lm and Da Neglected 
-200 + DX XCase 5: Da and Mw Neglected 
pnetedng Case 6: Lm, Da, and Mw Neglected 


Magnitude (dB) 


-270 


Phase (deg) 


-360 


-450 
1071 10° 101 10? 10° 104 10° 10° 
Frequency (rad/s) 


Figure 3. The bode plot of the EWB original model versus last 3 cases of ROMs 
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The accuracy of the proposed ROMs was further analysed using the time domain. For a simplified 
comparison, two groups of proposed models were made. Cases 1, 4, and 6 ROMs that showed different 
behaviour compared to the benchmark model are placed in the A group. Meanwhile, group B consists of the 
ROMs (2, 3, and 5) that offers a similar response to the benchmark. The clamping force is the output for the 
EWB model. The model output for each proposed model is plotted based on their group, according to the 
frequency-based accuracy result mentioned earlier. The clamping force response of each approximate model 
are as plotted in Figures 4(a) and (b). 

The step response shown in Figure 4(a) reveals that cases 1, 4, and 6 (group A) have a relatively 
large error. Two parts of the responses show errors from 0.1 to 0.3 seconds and 0.75 to 1.75 seconds, as 
highlighted in the graph. Meanwhile, the other group (group B) show very high accuracy responses which are 
close to the benchmark model, as in Figure 4(b). The model responses in the graph at 0.95 to 1.05 seconds 
showed no error spotted even though the graph was enlarged. 


Clamping Force Fn (N) 


EWB 5" Order Model 
— — — Case 2: Da Neglected 
= == =! Case 3: Mw Neglected 
muun Case 5: Da and Mw Neglected 


Clamping Force, Fn (N) 


EWB 5" Order Model 
— — — - Case 1: Lm Neglected 
=== =" Case 4: Lm and Da Neglected 
munun Case 6: Lm, Da and Mw Neglected 


0 0.5 1 1.5 2 2.5 3 3.5 0 0.5 1 1.5 2 2.5 3 3:5 
Time, t (s) Time, t (s) 
(a) (b) 


Figure 4. ROM clamping force step response for; (a) group A and (b) group B 


Based on the step response analysis, it is difficult to evaluate the accuracy of each model since the 
error is minimal and cannot be observed by normal observation. Therefore, the responses of both group 
models were investigated further by plotting the clamping force error, as in Figure 5. As displayed in 
Figure 5(a), the error for cases 1, 4, and 6 (group A) show the same magnitude and pattern. Figure 5(b) on the 
other hand, demonstrates the group B model deficient error, where it was found that the case 2 model 
accuracy is the best, followed by case 5 and case 3 models. The graph show that the error occurs only during 
the transient state and goes down to zero as time goes to infinity. 


x10 


Z Z 

É 2 
w w 
5 5 
È E 
w uw 
8 S 
6 © 
LW Ww 
D D 
£ £ 
Z i, 
E E 
S -350 S oo Seas ; 
O — — — - Case 1: Lm Neglected © 1/ Case 2: Da Neglected 

Abst ==" Case 4: Lm and Da Neglected -8 Y ==="! Case 3: Mw Neglected ] 
muuna Case 6: La, Da, and Mw Neglected aa Case 5: Da and Mw Neglected 
-450 L L L L L L L -10 
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Figure 5. ROM clamping force error for; (a) group A and (b) group B 
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In addition, the commonly used RMSE value was used to calculate the average error for the entire 
time range of the model response errors. The RMSE value and percentage for each approximated model are 
summarised in Table 2. The RMSE for case 2 was found to be the least with only 5.23e-12%, followed by 
case 5 and case 3 with 1.94e-5% and 2.80e-5%, respectively. Even though have a simpler system with almost 
zero error, the model reduction size produced in case 2 does not change. The best model reduction for case 5 
gives the simplest model with an acceptable error value. It is clearly shown that based on our design, only 3 
states are dominant which are motor position (8m), motor speed (©m), and motor current (Im). 


Table 2. The RMSE of clamping force for EWB reduce order ROMs 
Proposed model case number RMSE value (N) RMSE percentage (%) 


1 198.92 2.08 
2 5.00e-10 5.23e-12 
3 2.70e-3 2.80e-5 
4 198.92 2.08 
5 1.90e-3 1.94e-5 
6 198.92 2.08 


4. CONCLUSION 

The EWB parametric MOR is possible and can be made based on 3 main factors. The motor 
inductance, lead screw axial damping, and wedge mass are the factors that may have little influence on the 
overall system. The EWB widely used 5" full order model can be reduced up to 2™! order system, depending 
on the design of the EWB itself which has its parameters. Time and frequency domain response analysis is 
compulsory to determine which factor can be neglected and consequently choose the suitable ROM that can 
be applied. The RMSE can be used to select the best ROM. Even though the simplest EWB model is a 2™ 
order linear system, the best simplifications model may vary depending on EWB design parameters. 
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